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Abstract 

The concept of phase space amplitudes for systems with continuous degrees 
of freedom is generalized to finite-dimensional spin systems. Complex ampli- 
tudes are obtained on both a sphere and a finite lattice, in each case enabling a 
more fundamental description of pure spin states than that previously given by 
Wigner functions. In each case the Wigner function can be expressed as the star 
product of the amplitude and its conjugate, so providing a generalized Born in- 
terpretation of amplitudes that emphasizes their more fundamental status. The 
ordinary product of the amplitude and its conjugate produces a (generalized) 
spin Husimi function. The case of spin-i is treated in detail, and it is shown 
that phase space amplitudes on the sphere transform correctly as spinors under 
under rotations, despite their expression in terms of spherical harmonics. Spin 
amplitudes on a lattice are also found to transform as spinors. Applications are 
given to the phase space description of state superposition, and to the evolu- 
tion in phase space of the state of a spin-i magnetic dipole in a time-dependent 
magnetic field. 
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1 Introduction 



Weyl quantization (H is a one to one map between functions on a phase space T 
and operators acting on a Hilbert space T-L. Over the years the Weyl-Wigner (WW) 
transform W, which is the inverse of Weyl's quantization map, has been used to 
develop the phase space formulation of quantum mechanics that has successfully 
been applied to systems with continuous degrees of freedom I2H3]. 

There have been several attempts to generalize the Weyl correspondence to in- 
clude finite spin degrees of freedom within the WW formalism. The treatments 
of deGroot and Suttorp [4| and O'Connell and Wigner combined the continu- 
ous phase space picture with the finite spin degrees of freedom into a single WW 
transform. Buot [5| formulated a discrete WW transform on a periodic lattice phase 
space; Chumakov et al [6J defined a WW type quasi-probability function as a lin- 
ear combination of spherical harmonics on the phase space of the sphere S 2 ; and 
Berezin [7,8] introduced a general method of quantization by representing functions 
on S 2 in terms of covariant Q-symbols and contravariant P-symbols. 

The lattice phase space approach was further generalized to an array of N orthog- 
onal states to obtain a discrete Wigner function. For example Hannay and Berry [ 10 ] 
obtained a discrete Wigner function that resembled a delta-function on each of the 
(2N) 2 points of a 2N x 2N lattice phase space indexed by the set of integers modulo 
2N. Wootters |TT| generated a general class of discrete Wigner functions (applicable 
to spin systems) by constructing an N x N lattice phase space indexed by the set of 
integers modulo N, where N is a prime number, or a Cartesian product of lattices in 
the case of composite N. Gibbons et al [121 found a class of discrete Wigner functions 
on an N x N lattice phase space of a finite field composed of N k elements where N 
is prime and k £ Z + . Other approaches that dealt with spin systems but from a dif- 
ferent perspective are those of Chandler et al fll3fl who, motivated by the statistical 
characteristic function, generated three dimensional quasi-probability density func- 
tions for spin-i systems; and Leonhardt fl4 f , who applied precession tomography 
to spin systems to derive an expression for a discrete Wigner function. 

Phase space amplitudes have been defined for continuous degrees of freedom 
[T5ti2T| . They provide representations of quantum state vectors rather than density 
operators, and as such are more fundamental objects than Wigner functions M2"TH . 
In this paper we extend the concept of phase space amplitudes to finite-dimensional 
spin systems. We begin with a brief review of phase space amplitudes for continuous 
degrees of freedom and then modify the main results to define phase space ampli- 
tudes for a spin j system. Expressions for spin amplitudes and their correspond- 
ing Wigner functions are obtained on both the sphere S 2 and a lattice of dimension 
(2j + I) 2 . Our primary focus will be on spin-i amplitudes and explicit expressions 
for these will be given in terms of spherical harmonics and functions defined on 
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each of (2j + l) 2 = 4 lattice points. For convenience we consider all variables to be 
dimensionless and set Planck's constant H to unity. 



2 Phase space amplitudes 

In the phase space formulation of quantum mechanics, state vectors for a system 
with continuous degrees of freedom (p, q) only are represented by complex phase 
space amplitudes |[T5ti2T| . These are constructed from a pure state vector by 
introducing a fixed "window" state \ip) of unit length, to form the outer (dyadic) 
product 

* = |^)(y|- (1) 

For each choice of \tp), a set of distinct amplitudes for variable are then found as 
the images W( V I / ) under the WW transform, 



W(§)(p, q)/Vto = *(p, q) = Tr (§A(p, qfj /v 7 ^. (2) 



Here A(p, g) is the Weyl-Wigner-Stratonovich (WWS) operator kernel |22J. 

A phase space amplitude contains all essential information about a quantum 
state, and this representation of state vectors in phase space has many important 
properties. For instance, phase space amplitudes can be used to calculate expecta- 
tion values according to 

{A) = Tr(l$$ f ) — ► J V(p,q) [(A * *)(p, q)] dT, (dT = dpdq), (3) 

and transition amplitudes 

(H^) = f *i(p,q)Mp,q) dr, (4) 

where in particular 

<V#> = y *(p^)*(p, g) dr = l, (5) 

so that the quantity |^| 2 can be regarded as a quasi-probability distribution over T. It 
is remarkable |2~Tf that l^j 2 is the familiar Husimi distribution |24| if \tp) is chosen to 
be the ground state of the harmonic oscillator, and for other choices of the window 
state | tp), it is a generalized Husimi function (2TJ. In (O, the associative and non- 
commutative star-product appears, and in terms of phase space functions A(p,q) 
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and B(p, q), is given by |25] 



(A*B)(p,q) = (±\ jTr(A(p,q)A(j,\ q ')A(p'',q''))A(p',q')B(p'',q'') 

x dp'dq'dp"dq". (6) 

The Wigner function W(p, q) is defined for a pure state 1-0) as the image W(\ip)(ip\) 
and is commonly expressed as 11231 



9)/2tt 4— ► W(p, <l) = ^J e tpx i>(q + x/2)^{q - x/2) dx. (7) 

In analogy to the way that configuration space wave functions determine a prob- 
ability density, phase space amplitudes generalize the Born interpretation by deter- 
mining the Wigner function with the help of the star product 

W(p,q) = (y*y)(p,q), (8) 

which emphasizes the more fundamental status of these amplitudes. 



3 Spin amplitudes on the sphere S 2 

For each spin observable Jk, k = 1, 2, 3, represented by (2j + 1) x (2j + 1) her- 
mitian matrices with the fixed spin quantum number j £ 0, 1/2, 1 • • • , the inverse 
of Weyl's quantization map assigns a real valued function on the phase space of the 
sphere S 2 . However, pure spin states described in the Hilbert space %j = C 2jf+1 by 
eigenvectors \ j, m) of J3 and J 2 , 



m 



J, J 



3, 



(9) 



do not have images defined on5 2 . 

Within the framework of phase space amplitudes this problem is circumvented 
by first observing that spin states are vectors \tp) with 2j + 1 complex components 
i>m = (i) r ^IV') i n %j an d second, in a manner similar to the continuous case, by 
introducing a fixed spin state | (p) which allows one to extend the amplitude operator 
of ©to 



(10) 
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Normalization of | tp) requires that 



\<Pm\ 2 = l- (11) 



m=-j 



A basic feature of Weyl quantization is the existence of an operator valued func- 
tion or WWS kernel, the determination of which is essential for the definition not 
only of the Wigner functions on S 2 , as in the literature [6], but also of the spin am- 
plitudes in which we are interested here. A WWS kernel based on the axiomatic 
postulates of Stratonovich |22| and applicable to spin systems is known [26-28 J and 
given by 

#(0,0) = Zl m „\j,m")(j,m'\ (12a) 

m' ,m"=—j 

zi> m » = J^EE^-ir'C-ii^tw), d2b) 



2j + - 

1=0 m=-l 



where C^„ , r ] n are the standard Clebsch-Gordan coefficients, Y/ m are the spherical 
harmonics and ej are constants such that e J = 1 and ej = ±1. 

It is now a straightforward matter using ([10]) and ((121) , after making a suitable 
choice for the arbitrary fixed spin vector \ip), to define spin amplitudes on S 2 as 



*(0,<f>)=Tr [W(0,cf>) , (13) 



for each state \ tp) in ((TO)) . Whatever choice of fixed spin vector \tp) is made in ((10)) , the 
corresponding Wigner function on S 2 , which is defined as 

W(e,^ = Tr(\^){^(9,^), (14) 

is then expressed in terms of the spin amplitude by the star product 

W(M) = (¥**)(M). (15) 



There are various methods available to evaluate star products for spin systems [29 . 
|30| but for our purposes it will be appropriate to use the integral form by suitably 
modifying I© to get 

x m {9' , (t>')^>{6"^") dO'dO" , (16) 
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where dSl = sin OdOdcft is the invariant measure on S 2 . 

All functions defined on a sphere can be expanded in terms of a complete set of 
spherical harmonics, hence we can write the spin amplitude (TT3T > in the form 

2j I 

*(M)=EE ffl ^i). (it) 

Z=0 m=-l 

where, from the orthogonality property of the spherical harmonics, the coefficients 
are given by 

aim = [ [ *(M)*WM) sin0d0d<£. (18) 
Jo JO 

Formula ((17)) suggests that spin amplitudes will transform as tensors under rota- 
tions, but we will see that, to the contrary, these amplitudes transform as spinors 
when 2j is odd-integral, in particular when j = 1/2. 
A double application of ((17)) leads to 

|*(M)| 2 = E E ^mi^rna^WMWWM)- (19) 
h,mi h,m,2 

By analogy with the case of continuous variables |21], we can call 

|*(M)| 2 = <¥>|A'(M) QA j (e,ct>)\<p), (20) 

a generalized spin Husimi function if (9, <fi) on the sphere. In particular, if (up to an 
unimportant phase) | cp) is the Coherent Spin State (CSS) M3Tf 



m= -j V [J + m)\(j - m)l 

then (|20)) becomes the spin Husimi function [32 J. The product of two spherical har- 
monics with the same arguments can be written as a linear combination of single 
spherical harmonics in terms of the 3j-symbols, hence ((19)) can be re-expressed as 



im^ = EEEV (2 ' 1 + 1)(2 ': +1)(2 ' + 1) 



h,mi l2,m,2 l,m 



^ ~ a hmi a l2m2 



h h l\ h h l\ 7 r 

n o n ^ ( ^ (22) 

The Wigner function has a similar expansion, so that for any choice of fixed vector 
we can write 

W(9, 4>) = E W m (#, 0). (23) 
i=0 m=-Z 
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with coefficients 

r'ln f'TT 

km = 

JO JO 



W{9, 4>)Yi, m {e, <t>) sin 0d0#. (24) 



The symplectic group of transformations that arises in the case of continuous 
degrees of freedom, acting on functions of variables (p, q) in a 2f-dimensional phase 
space, is replaced by the group SU (2) of rotations acting on functions on the phase 
space of the sphere S 2 . Rotated spin amplitudes are generated from 

tt fl (M) = (W(£)*tf)(M), (25) 
where the rotation operator is given, for spin-i, by 

Rfr, a, (3) = e -i7^/2 e -ia? s /2 e -i^/2_ (26) 

For the case of spin-i systems (j = ^), setting the constants e l J 2 = +1, the WWS 
kernel ([121 reduces to 

A 1 / 2 (6', cf>) = (I + \/3 n ■ a)/2, (27) 
where / is a 2 x 2 unit matrix, a are the spin operators 

*-(::)•*-(: _0- 

and n = (sin # cos 4>, sin sin (j), cos 6*) is an arbitrary unit vector which is parameter- 
ized by the spherical polar coordinates 0<#<7r,0<(/>< 2tt. 

The spin-i amplitudes then follow by inserting ((10)) and (|27)l into (fl~3"|) to get 

l i 

= E tt ^(M), (29) 

i=0 m=-Z 

where the coefficients are given by 

«00 = ^1/2^1/2+^-1/2^-1/2, 

a w = V / 2V'-l/2^1/2> 

«10 = ^1/2^1/2-^-1/2^-1/2, 

an = --v/2^1/2^-1/2- ( 30 ) 

Using j22l , we find for the generalized spin Husimi function 

l i 

|tf (0, 0)| 2 = £ ft^m^- (31) 

i=0 m=-Z 
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with the associated coefficients 



boo = 1, 

Ki = -^(^-1/2^1/2 + ^-1/2^1/2)! 

ho = |*01/2| 2 |^1/2P - 1/2 1 2 1 V— X/2 1 2 . 

&11 = -^(^1/2^-1/2 + ^1/2^-1/2)- ( 32 ) 
The spin Husimi function follows from ||2"TT ) by choosing | cp) to be the CSS 

. . / cos hO \ .„_. 

I^)= iJ> ( 33 ) 

ye l<p sm^0J 

and observing that, in terms of expectation values of the spin operators, the elements 
of the density matrix are given by 

|V ± i/ 2 | 2 = U 1 ± ^1/2^1/2 = K(^) ± O 4 ) 

Thus, up on inserting the above into j20l ) leads to 11321 

|tf(0, ^)| 2 = H(0,<P) = const. (1 + ff (£)). (35) 

The corresponding Wigner functions are obtained by evaluating (fl6]> with ((27)) 
and ((29j), for any choice of fixed spin vector, to get 

1 1 

W(M) = V^£ E OrnYUW), (36) 

1=0 m=~l 

and its coefficients 



coo = 1, 

d_i = V / 2V'-l/2V'l/ 2 ! 

ClO = IV'l^l 2 - I V 7 — 1/2 1 2 , 

c u = -v / 2V ; i/2V'-i/2- (37) 

It is now possible to write down immediately the Wigner functions for any arbitrary 
spin vector by simply inserting appropriately chosen spin components. 

It has been established that spin amplitudes are functions on S 2 , and it only re- 
mains to be shown that they do indeed transform as spinors under rotations. To do 
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this consider a rotation through an angle a about the y-axis given by 

R = cos \al — i sin \oLG y . (38) 
The image W(R) of this rotation is given by the mapping 

W(cos \al — i sin ^aa y )(6, <j)) — > cos \a — iv^sm |asin#sin0. (39) 

Inserting this into ((25)) , and evaluating the star product, leads to an expression of the 
form 

l i 

* R (e,<t>) = V^Y,Y1 a lmYl, m (0,4>), (40) 
1=0 m =-l 

for the rotated spin amplitudes, with corresponding coefficients given by 

aoo = (^1/2^1/2 + ^-1/2^-1/2) cos \ a ~ (^-1/2^1/2 - ^1/2^-1/2) sin^a, 

an = ^(^-1/2^1/2 cosia + V1/2W1/2 Sirica), 

«io = (A/2<Pi/2 ~ ^-1/2^-1/2) cos|a- (^-1/2^1/2 + ^1/2^-1/2) sin 5«> 

an = -^(^1/2^-1/2 cos ia-V-1/2^-1/2 sin ^a). (41) 

Setting a = 27r in the above set of coefficients produces a change of sign in \P(0, 0). 
Therefore, even though the spin amplitudes are functions on 5 2 , expressed as linear 
combinations of spherical harmonics, they nevertheless transform as spinors. 

4 Spin amplitudes on the lattice 

Consider a quantum system characterized by two incompatible physical observ- 
ables A and B each with an associated finite number of (non-degenerate) eigenvalues 
a = {cti : i = —j, • • • , j} and /3 = {/% : i = —j, • • • , j} respectively. From these we 
construct a lattice consisting of (2j + l) 2 points (a, (3) in the plane. A generalized 
WW correspondence rule between observables Xj and their corresponding symbols 
X(a, (3) on each of the lattice points can be defined by 

X(a,/3) = Tv(x,-A(a,0)), (42b) 
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where A(a, j3) is a (2j + 1) x (2j + 1) lattice kernel matrix defined at the lattice points 
(a, /3) and required to satisfy the following properties 

A(a,/3) = A^a.jS) (43a) 

Tr(A(a,/3)) = 1 (43b) 

Tr(A(a,/3)A(c/,/3')) = (2j + 1) <W%' (43c) 

£A(a,/3) = (2j + 1) 1. (43d) 

a/3 

One possible form for the kernel A(a,/3) that guarantees the above properties are 
satisfied is H51ITTH 

A(a,P) = e~^ a '\a + a'/2)(a - a'/2\, a' = -j, ■ ■ ■ ,j, (44) 

a' 

where the |-)(-| are (2j + 1) x (2j + 1) matrices, whose elements are determined from 
the properties d43b - Alternatively, the matrix elements of A(a, j3) with respect to the 
basis are found from (1441 and given by the expression 

</3'|A(a,/3)|/3"> = e^rr^"-' 3 ') ^-^IM/J'+TOfl. (45) 

a' 

A double application of d42al) in conjunction with d42bD . leads to the star product 

(X*Y)(a, P) = | ^ X(a', /3')1>", 0")Tr (a(«, /3)A(a' ; /?') A(a", /?")) , 

(46) 

of two symbols defined on the lattice points. 

The image of the amplitude operator ^ of (flu) under the action of (|42b|l gives the 
amplitude on the lattice as 



*(a,/3)=TV(¥A(a,0)J, (47) 

with the corresponding Wigner function, obtained by evaluating the star product 
(|46]> with (dZll, given by 

W(a,/3) = — ^(tf**)^). ( 48 ) 
2 j + 1 

The lattice associated with spin j = \ consists of an array of four points desig- 
nated by {(0, 0), (0, 1), (1, 0), (1, 1)}. Defined on the lattice is a complex spinor am- 
plitude \& and an associated real Wigner function. The matrix elements for the lattice 
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kernel are in this case obtained directly from (|45)> and found to be 

The symplectic Fourier transform of the kernel A (a, (3), is 

D(a, p) = -L-Y l e^'-'^a', n (50) 

■I a'P' 

and in particular 

3(a, /3) = \ (A(0, 0) + (-l) a A(0, 1) + (-l)^A(l, 0) + (-l)^A(l, 1)) , (51) 
which leads to the Pauli operators 

3(0, 0) = I, 3(0, 1) = a x , 3(1, 1) = 5(1, 0) = a z . (52) 
Combining ((10)) and ((49)) to evaluate (|47)) we have for the spin amplitudes 

^(0,0) = ^1/2(^1/2 +^-1/2) + fC 1 + 0(^1/2^-1/2 -^-l/2^1/2) ; 

*(0,1) = ^-1/2(^-1/2 + W1/2) + 2K 1 - 0(^1/2^-1/2 -^-1/2^1/2)1 

#(1,0) = Tpi/ 2 (l/>l/2 ~ ^-1/2) - |(1 + 0(^1/2^-1/2 - ^-1/2^1/2)1 

#(1,1) = '0-1/2(^-1/2 -^1/2) - |(1 - 0(^1/2^-1/2 -^-1/2^1/2)- ( 53 ) 

These lead immediately to 

|#(0, 0)| 2 = 1(^/2 + ^-1/2)^1/2 + 2(^1/2 + ^-1/2)^1/2 1 2 /2, 

l#(0, 1)| 2 = 1(^1/2 + ^-l/ 2 )?-l/2+ »(^l/2 +^-l/2)^l/2| 2 /2, 
|#(1, 0)| 2 = -^-1/2)^1/2 + »(^l/2 -^-1/2)^1/2 iV 2 , 

l#(l, !)! 2 = 1(^1/2 -^-1/2)^-1/2 + ^1/2 -^-i/2)Vi/2| 2 /2, (54) 

which defines a discrete spin Husimi distribution |33| on the lattice. Using the am- 
plitudes ((53)l in (|48)) it is found that, whatever choice is made for the fixed spin state 
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\tp), the Wigner function for any given state \ip) is 

W(0,Q) = [?l /2 (^l/2 + i>-l/2) + + 0(^1/2^-1/3 - ^-l/2?l/ 2 )]/2, 
W(0, 1) = [^-1/2^-1/2 + V^l/2) + 5(1 " »)(^l/2?-l/2 " V'-l/2?l/ 2 )]/2, 
W(l, 0) = ^ 1/2 (^l/2 " ^-1/2) " l(l + 0(^1/2^-1/2 " ^-l/2?l/ 2 )]/2, 

W(l,l) = [V--i/2(?-i/ 2 -?i/ 2 )-K 1 - i )M-i/ 2 -^i/2?i/ 2 )]/2. (55) 

Note that even on a lattice the Wigner function can take negative values and there- 
fore cannot be interpreted as a probability density function. This can be seen by 
choosing, for example, ^1/2 = (1 — O/v 7 ^, ip-1/2 = V2(l + 0- 

The images W(R) of a rotation about the y-axis, given by l(38)l , on each of the four 
lattice points are then given by the relations 

#(0,0) =12(1,1) = e~ ia/2 , R(0,l)=R(hO) = e ia/2 . (56) 

Using (|56)) and ((53)) , the rotated spin amplitudes on each of the lattice points are 
obtained from the formula 

9 R (a,P) = (#* *)(«,/?). (57) 
This is evaluated according to the formula ((46j) and yields 

*r(0,0) = ^1/2(^1/2 + ^-1/2) cos|a + ^1/2(^1/2 -^-1/2) sin ia + |(l + 

x ((^1/2^-1/2 - ^-1/2^1/2) cos 3 a - (^1/2^1/2 + ^-1/2^-1/2) sin 3«)> 
$r(0,1) = ^-1/2(^-1/2 + ^1/2) cos ±a + ^1/2(^1/2+ W-1/2) sinia + i(l-i) 

x ((^1/2^-1/2 - ^-1/2^1/2) cos 2« - (^1/2^1/2 + ^-1/2^-1/2) sin 2 a )' 
$r(1,0) = ^1/2(^1/2 - ^-1/2) cos 5 a -^1/2(^1/2 + ^-1/2) sin |a - |(1 + 

x ((^1/2^-1/2 - ^-1/2^1/2) cos |a - (^1/2^1/2 + ^-1/2^-1/2) sin §a), 
*r(1,1) = ^-1/2(^-1/2 -^1/2) cos §a + ^1/2(^-1/2 -^1/2) sin ~a - ~(1 - i) 

x ((V'l/2^-1/2 - ^-1/2^1/2) cos §a - (^1/2^1/2 + ^-1/2^-1/2) sin 

(58) 

Setting a = 2-7T in the above set of amplitudes results in a change of sign in *&(a, /3) 
and thus shows that spin amplitudes defined on the lattice also transform as spinors 
under rotations. 
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5 State Superposition in Phase Space 



The most important new feature that amplitudes bring to the description of quan- 
tum systems in phase space is their linear superposability, a property that they in- 
herit from the superposability of quantum state vectors by virtue of the definitions 
(f!3l > and d47b , which are linear in the relevant state vectors. Because the definition 
((151) and ((48|) of a Wigner function is not linear in the corresponding state vector, 
there is no simple way to construct W from W\ and W2 using phase space methods, 
when W, W\ and W2 are Wigner functions corresponding to the pure states \ip), \ip\) 
and \ip2) in the superposition 

|V> =a|Vi>+/3|V 2 ). (59) 

However, if the phase space amplitudes ^1 and *2 corresponding to |?/>i) and \ip2) 
are constructed, for one particular choice of reference (window) state vector \ <p), then 
we evidently have for the amplitude \P corresponding to the superposition j59l , sim- 
ply 

q> = aV 1 + /3* 2 . (60) 

Then the corresponding Wigner function W can be constructed from the generalized 
Born interpretation, as 

W = ^ *^ = (a*! + /3* 2 ) * (a*i + ^2) 

= |a| 2 *i + |/3| 2 *2 + a/3*i **2 + /3a^ 2 **i 

= |a| 2 Wi + \/3\ 2 W 2 + a/?*i**2~ + /3a* 2 **i- (61) 

The cross terms in the final expression here emphasize the importance of spin am- 
plitudes, over and above that of Wigner functions, in the phase space description of 
superposition states. 

These considerations apply to amplitudes defined on the sphere or the lattice, just 
as they do in the continuous case |2TH . Similar remarks apply also to the construction 
of spin Husimi functions: if H, Hi and H2 correspond to \ip), \ipi) and ^2) in (|59)l , 
then there is no straightforward phase space method to construct H from H\ and H2. 
However, with the introduction of phase space amplitudes, on the sphere or lattice, 
we have 

H=\$\ 2 = \a^ 1 +f3^ 2 \ 2 

= |a| 2 #i + \(3\ 2 H 2 + a/3#i*2" + 0a¥ 2 *i- (62) 

Note that * carries the same information about the relative phases of \ipi) and 1^2) as 
does the linear combination \ tp). Interference phenomena will therefore be described 
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transparently in phase space in terms of amplitudes, whereas the description of such 
phenomena in terms of Wigner or Husimi functions is necessarily less clear. 

6 Example: Spinor Amplitudes in a Magnetic Field 

An important example of a quantum control system involves the manipulation 
of a spin-i magnetic dipole with an external magnetic field. Procedures of this kind 
are prevalent in nuclear magnetic resonance (NMR) spectroscopy. In this case non- 
interacting spin-i dipoles are immersed in an external homogeneous magnetic field 
and perturbed from equilibrium by a controlled magnetic field oscillating in the 
plane perpendicular to the external field. It is illuminating to investigate the form 
and behaviour of a spinor amplitude representing an isolated spin-i dipole under 
the influence of an NMR control environment. 

Consider then such a dipole exposed to an external constant magnetic field B z 
acting along the z-axis and subjected to a (weaker) controlled magnetic field B xy 
rotating in the xy-axis. The governing time dependent Hamiltonian at the time t 
when a phased sinusoidal pulse is applied is given by |34| 

H = \ujoo z + |cj nu t [cos(w rcf t + Xp)% + sin(w ref t + Xp)°y] , (63) 

where uq = —jB z is the Larmor frequency uj nut = —jB xy is the nutation frequency 
w rc f is the spectrometer reference frequency 7 the gyromagnetic ratio, \p is the- phase 
angle, and a x ,a y , d z are the Pauli matrices as in ((28]). 

A simplification is obtained by eliminating the time from the Hamiltonian <(63j) by 
transforming to a frame which is rotating at the Larmor frequency. Under this trans- 
formation, the time independent Hamiltonian at the time of the sinusoidal pulse 
takes the form 

H rot = \uj ms a z + 5 w nut (a x cos x P + Vy sin x P ), (64) 

where oj res = ujq — w re f is the resonance off -set. The spinor amplitude derived in this 
frame will appear to be stationary. 

The amplitude operator (fTOb , with j = \, in the rotating frame is found from 

%' = f7#, (65) 

where the evolution operator U for this Hamiltonian is 

U = e~ iSmtt , (66) 
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and is expressed in matrix form as 



- _ /cos \a - iu res sin ^a/uj cS -ie iXp w nut sin \a/u eS 
V —ie lx p uj nut sin ^ a/ w e ff cos | a + iw res sin | a/ w e g 

Here the flip angle q = u; c ff i and the quantity w c ff = v ^rcs + w nut * s the magnitude 
of the off-resonance rotation frequency that gyrates around the tilted axis. 

The spinor amplitude ^(9,(j)) in the rotating frame can be found by evaluating 

(U*V)(0,<j>) or Tr(U$A 1/2 (6, <j>)). (68) 

Expressed in terms of spherical harmonics the spinor amplitude has the form (|29)) , 
with coefficients now given from d68l by 

«00 = Ol/2^1/2 +^-1/2^-1/2) COS \ a ~ ^[(^1/2^1/2 - ^-l/2^-l/2) w rcs 

+ (e" lXp V ; -i/2^i/2 + e lXp V'i/2^-i/2) w iiut] sin ±a/w cfr , 

= V / 2(^_i/ 2 ^l/2 cos ^« + i[V ; -l/2^1/2 w res 

- e iXj, Vi/2^i/2 w nut]sinia/u; cfr ), 

«10 = (^1/2^1/2 - ^-1/2^-1/2) COS 2"° ~ ^[(^1/2^1/2 +^-1/2^-1/2)^ 

+ (e _% ^-i/2^i/2 - e lXj, Vi/2^-i/2) w nut] sin ia/w efT , 
an = -^(^1/2^-1/2 cos ia - i[^i/2^-i/2^rcs 

+ e" JXp V-i/2^-i/2^nut] sin ia/wgfif). (69) 

As resonance is approached, that is when wo ~ <^rcf/ we find a = w nut t and 
w e ff = w nut and the coefficients l|69)l reduce to 

«oo = (^1/2^1/2 + ^-1/2^-1/2) cos|a 

-i(e" lXp V-i/2^i/2 + e lXp ^i/ 2 W-i/2) sin 5a, 

Ow = ^2(^-1/2^1/2 cos|ci!-'je % ^i/2^i/2Sin|a), 

«io = (^1/2^1/2 -^-1/2^-1/2) cc, s \ol 

-i(e~" lXp if>-i/2<Pi/2 ~ e tXp Vi/2^-i/2) sin 

an = -V2(il)i/2<P-i/2 c °s 5« - se~ lX3 "0-i/2^-i/2 sin 5")- (70) 

Thus, at resonance the spinor amplitude is subjected to a phased rotation through an 
angle a, where the amount of rotation is dictated by the duration of the controlled 
magnetic field. In particular, when the phase Xp = then ( |70l are reduced to the 
coefficients that represent a rotation about the y-axis through an angle a (compare 
with d4Tb). By switching off the magnetic field, that is setting w nut = in ( fTDl , the 
spinor amplitude returns back to its initial equilibrium position of precessing about 
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the z-axis. 

In the case of a lattice with an associated spin j = |, the phase space is com- 
prised of an array of four points. The lattice spinor amplitude for an NMR type con- 
trol configuration for a single spin-^ system that is described by the rotating frame 
Hamiltonian d64t is obtained by using d65l ) in j47l ) to give 

*(a,P)=Tr(u*A(a,pj), (71) 

where U, as before, is given by l|67)l , and the lattice kernel A by ((49} . Evaluating ((71) , 
we find for the lattice spinor amplitude 

#(0,0) = ^1/2(^1/2 + ^-1/2) cos \a 

+ (^-1/2 - ^1/2) (*0i/2 w res + ^-i/2^nute" lXp ) sin ^a/a; cff 

+ l(! + O((01/2^-l/2 - ^-1/2^1/2) cos 5 Q 
+ (00-1/2^1/2 - 01/2^-1/2)^08 

- (^1/2^1/2 e * Xp +^-i/ 2 ^-i/2 e ^ Xp ) w nut)sinia/w cff ), 
#(0,1) = ^-1/2(^1/2 + ^-1/2) cos 

- (<?l/2 - W-l/2) (V , -l/2^res - ^l^mite 1 **) sin ia/w ofr 
+ 5(1 + 0(00-1/2^1/2 - ^1/2^-1/2) cos §a 

+ (00-1/2^1/2 - 01/2^-l/2Vres 

- (^1/2^1/2 e * Xp + V , -i/2^-i/2 e ~ 4Xp )^nut)sinia/a; eff ), 
#(1,0) = ^1/2(^1/2 -1P-1/2) cos |a 

- (^-1/2 + ^1/2) (*0i/2^res + 0-i/2^nute~ iXp ) sin ^a/a; cff 

- 5(1 + ((^1/2^-1/2 - ^-1/2^1/2) cos |a 

- (001/2^-1/2 - 0-1/2^1/2)^8 

+ (^1/2^1/2 e * Xp + V'-i/2^-i/2 e ^ Xp ) w nut)sinia/w eff ), 
#(1,1) = ^-1/2(^-1/2 -^1/2) cos 5a 

+ (^1/2 + W-1/2) (^-i/2 w res - 0i/ 2 Wnute iXp ) sin ^a/w ofr 
+ 1(1 + ^((^1/2^-1/2 - ^-1/2^1/2) cos \a 

+ (001/2^-1/2 - 0-l/2^1/2) w res 

+ (^1/2^1/2 e ' XP + ^-l/2^-l/2 e ~* XP )^nut) Sm ±a/u eS ), (72) 

where the flip angle a = oj e gt and oj e s = \/w 2 es + w„ ut . When u;q ~ u; re f, near reso- 
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nance, we find that ((72|) , after a slight rearrangement, reduce to 



#(0,0) = ^1/2(^1/2 +^-1/2) cos la + tp_ l/2 e'' lXp (Tp_ l/2 - up 1/2 ) sin±a 

+ ki 1 + 0((^l/2^-l/2 - ^-1/2^1/2) COS |q 

- (^1/2^1/2 e % +li/2^-i/ 2 e"' & )sm^), 

#(0, 1) = ^-1/2(^-1/2 + ¥1/2) cos \a + ^-1/2(^-1/26"^ + iVi/2 e% sin | a 

+ sC 1 _ ((^1/2^-1/2 - ^-1/2^1/2) cos |a 

+ (^1/2^1/2 e* Xp +V'-l/2^-i/2e - * Xp )sin^a), 
#(1,0) = ^1/2(^1/2 -V'-i/2)cosia-V-i/2e" JXp (^_i/2 + ^i/2) s i n ia 

- 5(1 + ((^1/2^-1/2 - ^-1/2^1/2) cos \a 

- (^1/2^1/2 e tXp +^_i/ 2 ^_i/2e~ 4Xp )smia), 

#(1,1) = ^-1/2(^-1/2 -^i/2)cosia + ^_ 1/2 ('0-i/2e" lXp - #i/ 2 e iXp )sinia 

- 5 (! - 0((^l/2^-i/2 - ^-1/2^1/2) cos ia 

+ (V»i/2^i/2e <Xp + ^i/2^-i/ 2 e~ 4Xp ) sin ±a). (73) 

At resonance the lattice spinor amplitude experiences a phased rotation through an 
angle a, which is determined by the length of time that the magnetic field is applied. 
Choosing a phase angle of Xp = vr/2 in (f73]> results in a spinor amplitude repre- 
senting a rotation about the y-axis through an angle a as in (|58)) . By turning off the 
magnetic field, that is by setting oj nut = 0, the lattice spinor amplitude returns to its 
equilibrium state. 

7 Concluding remarks 

We have extended the concept of phase space amplitudes to finite spin systems 
and introduced the notion of a spin amplitude to represent pure spin states in phase 
space. Spin amplitudes and their associated Wigner functions have been expressed 
as linear combinations of spherical harmonics on the sphere. In addition, a general- 
ized Weyl correspondence has been adapted to define spin amplitudes and Wigner 
functions on a lattice phase space. We have concentrated on amplitudes for a single 
fixed spin, in particular a spin-i system. Our approach can be extended in the ob- 
vious way to define amplitudes for multipartite spin systems on Cartesian products 
of one-particle phase spaces, either spheres or lattices. 

As an example the general theory has been applied in detail to the case of spin- 
\, revealing some important features of the approach. The liberty to choose any 
normalized fixed spin reference state \ip) introduces a degree of arbitrariness into 
the definition of the spin amplitudes, analogous to the freedom to choose a "window 
state" in the continuous case [21|. It is an important problem to optimize this choice 
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for a given quantum spin system. We emphasize however, as can be seen from ((HJ) 
and ((48|) , that the associated Wigner functions are independent of this choice, just as 
in the continuous case. 

The spin amplitudes are completely described in terms of a combination of sys- 
tem and fixed state spinor components and these amplitudes transform as spinors 
under rotations on both the sphere and lattice thus further supporting their funda- 
mental status. Spin amplitudes in phase space can be superposed like state vectors, 
unlike Wigner functions or Husimi functions. They enable an expanded phase space 
description of quantum spin systems, as in our example of a magnetic dipole in a 
time-dependent magnetic field. 

It is our view that the representation of spin states in phase space holds the 
promise of new physical insights and a novel perspective in areas like quantum com- 
puting, quantum control systems, and quantum information. We hope to return to 
some of these aspects in the future. 
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